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Abstract 

A simple and transparent example of a non-autonomous flow system, with hyper- 
bolic strange attractor is suggested. The system is constructed on a basis of two coupled 
van der Pol oscillators, the characteristic frequencies differ twice, and the parameters 
controlling generation in both oscillators undergo a slow periodic counter-phase varia- 
tion in time. In terms of stroboscopic Poincare section, the respective four-dimensional 
mapping has a hyperbolic strange attractor of Smale - Williams type. Qualitative rea- 
soning and quantitative data of numerical computations are presented and discussed, 
e.g. Lyapunov exponents and their parameter dependencies. A special test for hy- 
perbolicity based on statistical analysis of distributions of angles between stable and 
unstable subspaces of a chaotic trajectory has been performed. Perspectives of further 
comparative studies of hyperbolic and non-hyperbolic chaotic dynamics in physical 
aspect are outlined. 

Mathematical theory of chaotic dynamics based on rigorous axiomatic foundation ex- 
ploits a notion of hyperbolicity, which implies that all relevant trajectories in phase space 
of a dynamical system are of saddle type, with well defined stable and unstable directions 
[H 12 HH HI- Hyperbolic systems of dissipative type, contracting the phase space volume, 
manifest robust strange attractors with strong chaotic properties. The robustness (structural 
stability) implies insensitivity of the motions in respect to variations of equations governing 
the dynamics. In particular, positive Lyapunov exponent depends on parameters in smooth 
manner, without flops into negative region characteristic to non-hyperbolic attractors. A 
Cantor-like structure of the strange attractor persists without qualitative changes (bifurca- 
tions), at least while the variations are not too large. Textbook examples of these robust 
strange attractors are represented only by artificial mathematical constructions associated 
with discrete-time models, e.g. Plykin attractor and Smale - Williams attractor (solenoid). 

It seems that the mathematical theory of hyperbolic chaos has been never applied con- 
clusively to any physical object, although concepts of this theory are widely used for in- 
terpretation of chaotic behavior of realistic nonlinear systems. On the other hand, feasible 
nonlinear systems with complex dynamics, such as Lorenz and Rossler equations, chaotic 
self-oscillators, driven nonlinear oscillators etc. do not relate to the true hyperbolic class 
|H El E]- As a rule, observable chaos in these systems is linked with a so-called quasiat- 
tractor, a set in phase space, on which chaotic trajectories coexist with stable orbits of 



high periods (usually, they are non-distinguishable in computations at reasonable accura- 
cies). Mathematical description of quasiattractors remains a challenging problem, although 
in physical systems the non-hyperbolicity is masked effectively due to presence of inevitable 
noise. In few cases, e.g. in Lorenz model in some appropriate domain of parameter space, 
dynamics is proved to be quasi-hyperbolic (with some restrictions concerning violation of 
smoothness conditions) [7j. 

I am aware of only two theoretical works, which discuss examples of true hyperbolic 
dynamics in context of systems governed by differential equations. One relates to a system 
called triple linkage, which allows in a frictionless case a description in terms of orbits on 
a surface of negative curvature. In dissipative case, it gives rise to a hyperbolic chaotic 
attractor 8J. Another work deals with a 3D flow system motivated by neural dynamics and 
argues in favor of existence of an attractor of Plykin type in the Poincare map associated 
with the flow p. 

In this Letter, I suggest an essentially simpler and transparent example of a non-autonomous 
flow system, which apparently manifests a hyperbolic attractor. In terms of stroboscopic 
Poincare map, it is an attractor of the same kind as the Smale - Williams solenoid, but 
embedded in a 4D rather then 3D state space. 

The system is constructed on a basis of two van der Pol oscillators with characteristic 
frequencies uq and 2ujq, respectively. The control parameters of the oscillators responsible 
for the Andronov - Hopf bifurcations in the autonomous subsystems are forced to swing 
slowly, periodically in time. On a half-period, the first oscillator is above the generation 
threshold, while the second one is below the threshold. On another half-period, a situation 
is opposite. Next, we assume that the first oscillator acts on the partner via a quadratic term 
in the equation. The generated second harmonic component serves as a primer for the second 
oscillator, as it comes off the under-threshold state. In turn, the second oscillator acts on 
the first one via a term represented by a product of the dynamical variable and an auxiliary 
signal of frequency ujq. Thus, a component with the difference frequency appears, which fits 
resonance range for the first oscillator and serves as a primer as it starts to generate. 

Summarizing this description, we write down the following equations: 



where x and y are dynamical variables of the first and the second oscillators, respectively, 
A is a constant designating amplitude of the slow swing of the control parameters, e is a 
coupling parameter. 

We assume that the period of swing T contains an integer number of periods of the 
auxiliary signal: T = 2ttN/ujo. Thus, our set of non-autonomous equations has periodic 
rather than quasiperiodic coefficients. It is appropriate to treat the dynamics in terms of 
stroboscopic Poincare section using a period-T sequence of time instants. The Poincare map 
is four-dimensional and acts in a space of vectors {x, x/uo, y, y/(2uio)}- 

The system (JTJ) operates as follows. Let the first oscillator have some phase (p on a 
stage of generation: x cx cos(co> £ + V 9 )- Squared value x 2 contains the second harmonic: 
cos(2cl> £ + 2<p), and its phase is 2<p. As the half-period comes to the end, and the second 
oscillator becomes excited, the induced oscillations of the variable y get the same phase 2<p. 
Mixture of these oscillations with the auxiliary signal transfers the doubled phase into the 
original frequency range. Hence, on the next stage of excitation the first oscillator accepts 
this phase 2ip too. Obviously, on subsequent stages of swing the phases of the first oscillator 
follow approximately the mapping 



x — (Asm2irt/T — x 2 )x + uj^x = sysmuiot, 
y - (-A sin 2-Kt/T - y 2 )y + Atu^y = ex 2 , 




y? n+ i = 2ip n (mod27r). 



(2) 



Figure 1 shows a typical pattern of time dependence of x and y from numerical solution 
of Eqs. (0) by Runge - Kutta method for particular parameter values ujq = 2ir, T = N = 
10, A = 3, e = 0.5 together with a diagram of empirical mapping for phase (p n+ \ versus <p n . 
The phases are determined at the centers of the excitation stages for the first oscillator: 

_ J arctan (lu^x/x), x > 0, , . 

1/9 1 arctan (u^x/x) + 7r, x > 0, 

and are plotted over a sufficiently large number of the basic periods T. The mapping for the 
phase looks, as expected, topologically equivalent to the relation (J2J). (Some distortions arise 
due to imperfection of the above qualitative considerations and of the definition of phase; the 
correspondence becomes better at larger period ratios N.) Chaotic nature of the dynamics 
reveals itself in a random walk of humps in respect to the envelope of the generated signal 
on subsequent periods of swing. 




Figure 1: A typical pattern of time dependence for variables x and y obtained from numerical 
solution of Eqs. (JJJ for u = 2n, T = N = 10, A = 3, e — 0.5 (a) and a diagram of empirical 
mapping for phases of the first oscillator defined in the centers of the stages of excitation 
numbered by n 

In terms of stroboscopic Poincare map, attractor of the system corresponds exactly to 
the construction of Smale and Williams. In the four- dimensional state space, the direction 
associated with the phase (p is expanding and gives rise to Lyapunov exponent estimated 
as Ai ~ T _1 log 2. Three rest directions are contracting, and they correspond to a three- 
dimensional stable manifold of the attractor. Three respective Lyapunov exponents are 
negative. Interpreting the stroboscopic Poincare mapping, we may imagine a solid torus 
embedded in 4-dimensional space and associate one iteration of the map with longitudinal 
stretch of the torus, with contraction in the transversal directions, and insertion of the doubly 
folded "tube" into the original torus interior. 

In computations, the Lyapunov exponents were evaluated with a help of Benettin's al- 
gorithm pm ITT] from simultaneous solution of Eqs. (0) together with a collection of four 
exemplars of the linearized equations for perturbations: 

x + 2xxx — (Asm27it/T — x 2 )x + u^x = eysmuj Q t, , , 

y + 2yyy - {-A sin 2nt/T - y 2 )'y + 4u%y = 2exx. 

In a course of the solution, at each step of the integration schema, the Gram - Schmidt or- 
thogonalization and normalization were performed for four vectors jx(t), x(t)/uo, y(t),y(t)/2(jQ 
and the mean rates of growth or decrease of the accumulated sums of logarithms of the norms 



(after the orthogonalization but before the normalization) were estimated. As found, the 
Lyapunov exponents for the attractor at the above mentioned parameters are A x w 0.068 ~ 
T- 1 log2, A 2 « -0.35, A 3 « -0.59, A 4 « -0.81. 

If the attractor is indeed hyperbolic, the chaotic dynamics must be robust and retain its 
character under (at least small) variations of the equations. As checked, this is indeed the 
case. In particular, the largest Lyapunov exponent is almost independent on parameters, 
and the rest of them manifest regular parameter dependences, as seen in Fig. 2. The left 
edge of the diagram corresponds to violation of the hyperbolicity. 

Dynamical behavior of the same kind is observed at other integer period ratios N, includ- 
ing essentially smaller ones, e.g. N=4. Figure 3 shows portrait of the strange attractor in the 
Poincare section in projection onto the plane (x, x) at u = 2ir, T = N = 4, A = 8, e — 0.5. 
It looks precisely as the Smale - Williams attractor should look like. Observe fractal transver- 
sal structure of "strips" constituting the attractor. For this attractor the Lyapunov expo- 
nents are Ai « 0.167, A 2 ~ —0.72, A 3 « —1.03, A 4 « —1.50. An estimate for fractal 
dimension from Kaplan - Yorke formula yields D w 1.23, and that from the Grassberger - 
Procaccia algorithm is D w 1.26. 




Figure 2: Computed Lyapunov exponents of the system versus parameter A at ujq = 2tt, 
N = T = 10, e = 0.5. Observe that the largest exponent remains almost constant in the 
whole interval of hyperbolicity being in good agreement with the estimate Ai = T" 1 log 2. 



It is interesting to perform a direct numerical test for hyperbolicity of the attractor. Idea 
of such test was suggested in Refs. [T2| and [T3] and applied for verification of hyperbolicity 
of trajectories of dynamical systems, which have one stable and one unstable directions. The 
procedure consists in computation of vectors of small perturbations along the trajectory in 
forward and inverse time with measuring angles between the forward-time and backward- 
time vectors at points of the trajectory. If zero values of the angle do not occur, i.e. the 
statistical distribution of the angles is essentially separated from zero, one concludes that the 
dynamics is hyperbolic. If the statistical distribution shows non- vanishing probability for zero 
angle, it implies non-hyperbolic behavior because of presence of the homoclinic tangencies 
of stable and unstable manifolds. In dissipative case these tangencies are responsible for the 
occurrence of quasiattractor. 




Figure 3: Portrait of the strange attractor in the stroboscopic Poincare section t — 1 (mod 
N) in projection onto the plane (x, x) at uj® = 2n, T = N = 4, A — 8, e — 0.5 

In our system only unstable subspace is one-dimensional, and the stable one is 
three-dimensional. Therefore, the method needs a modification. An adopted algorithm 
consists in the following. First, we generate a sufficiently long representative orbit on the 
attractor {x(t), x(t)/uj , y(t), y{t)/2oJo\ from the numerical solution of Eqs. (0). Then, we 
solve numerically the equations for a perturbation forward in time. In a course of 
the solution, normalization of the vector a(t) = [x(t),x(t)/u ,y(t),y(t) /2u } is performed 
after each step of integration to exclude the divergence. This vector determines an unstable 
direction at each point of the orbit. Next, we solve a collection of three exemplars of equations 
(JIJ) in backward time along the same trajectory {x(t), x(t)/u>o, y(t), y(t)/2ujo} to get three 
vectors {b(t), c(t), d(t)}. To avoid dominance of one eigenvector and divergence, we use 
the Gram - Schmidt orthogonalization and normalization of the vectors at each step of the 



numerical integration. Now, at each point of the trajectory, all possible linear combinations 
of {b(i), c(t), d(t)} define a three-dimensional stable subspace of perturbation vectors. 

To estimate an angle a between the one-dimensional unstable subspace and the three- 
dimensional stable subspace we first construct a vector v(t) orthogonal to the three-dimensional 
subspace, with components determined from a set of linear equations v(t) ■ h(t) = 0, 
v(t) ■ c(t) = 0, v(t) • d(t) = 0. Then, we compute an angle ft G [0,7r/2] between the 
vectors v(t) and a(£): cos ft = |v(t) • c(i)|/|v(£)||c(£)| and set a = ft — n/2. 




Figure 4: Histograms for distributions of angles a between the stable and unstable subspaces 
for the system (0) with u = 2tc, e = 0.5 obtained from computational procedure described 
in the text: (a)JV = 10, A = 3 and (b) JV = 4, A = 8 

Figure 4 shows histogram for the distribution of angles a between the stable and unstable 
subspaces for the system (JTJ obtained from computations at the two mentioned sets of 
parameter values. Observe clearly visible separation of the distributions from zero values of 
a. So, the test confirms hyperbolicity of the attractors. 

In spite of simplicity of the presented example, I believe it is significant as a feasible 
system, which may be designed as a physical device, e.g. on a basis of two interacting 
electronic oscillators. It opens an opportunity for experimental studies of hyperbolic chaos 
and its features predicted by the mathematical hyperbolic theory (robustness, continuity of 
the invariant measure, insensitivity of statistical characteristics of the motions in respect 
to noise, etc.). In addition, it makes conclusive comparative examination of dynamics of 
hyperbolic and non-hyperbolic systems. 

In a sense, breakthrough into the hyperbolic domain is a decisive step. Now one can 
construct many other examples of systems with hyperbolic attractors: Because of robustness 
of such attractor, any variation of the right-hand parts of the equations will not destroy the 
hyperbolicity, at least while they are not too large. Apparently, in this way it is possible 
to design examples of autonomous systems with hyperbolic strange attractorsb as well, via 
modification of the system supplimenting additional equations for dynamical variables, which 
would represent the swing and the auxiliary signals. 
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